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In  recent  years  gadolinium  scandium  gallium  garnet,  GdjScjGajO^ 
(GSGG),  has  received  growing  attention  as  a  laser  host  material  1 1-1 1 1. 
The  luminescence  and  stimulated  emission  from  trivalent  neodymium  in 
GSGG  were  first  reported  by  Kaminskii  and  coworkers  1 1].  Struve  et  al 
[2]  observed  broad  fluorescence  from  GSGG  doped  with  trivalent 
chromium.  They  were  able  to  demonstrate  Cr3+  laser  action  in  their  crys¬ 
tals.  Several  groups  [3-7 j  discovered  that  co-doping  GSGG  with  Nd3+ 
and  Cr3+  resulted  in  more  efficient  energy  transfer  to  the  Nd3+  energy  lev¬ 
els  involved  in  laser  action.  The  spectroscopic,  optical,  and  ther¬ 
momechanical  properties  of  Nd3+:Cr3+:GSGG  recently  reported  by 
Krupke  and  coworkers  [  1 1 J  suggest  potential  for  this  material  as  a  power¬ 
ful  laser. 

Several  years  ago  we  began  a  theoretical  analysis  of  the  crystalline  electric 
field  (CEF)  splitting  of  the  energy  levels  of  trivalent  neodymium, 
gadolinium,  and  chromium  in  GSGG  [12-14).  At  that  time  only  the  ex¬ 
perimental  CEF  splitting  of  the  [4/|7  and  [4F]3/2  manifolds  of  trivalent 
neodymium  had  been  reported  [1,15,16].  Recently  Struve  and  Huber  [  17| 
published  their  work  on  the  energy  levels,  line  broadening,  and  transition 
probabilities  based  on  the  optical  spectra  of  Cr3+:GSGG.  Using  as  an  ap¬ 
proximation  the  cubic  energy  level  scheme  of  Sugano  and  Tanabe  [  18, 19|, 

A  /J 

they  concluded  that  the  T2  level  is  admixed  into  the  E  level  via  spin- 
orbit  coupling  [20 j.  However,  no  CEF-splitting  calculations  addressing 
the  correct  symmetry  for  Cr3+,  Nd3+,  or  Gd3+  in  GSGG  have  been 
reported. 

This  paper  presents  the  absorption  spectra  and  energy  levels  for  Gd3+, 
Nd3+,  and  Cr3+  in  GSGG  up  to  35,000  cm  1  recorded  at  liquid  helium,  liq¬ 
uid  nitrogen,  and  room  temperatures.  The  CEF-splitting  calculations  for 
Gd3+,  Nd3+,  and  Cr3+  are  based  on  the  complete  diagonalization  of  a 
Hamiltonian  including  Coulombic,  spin-orbit,  and  CEF  terms  in  symmetry 
appropriate  to  the  trivalent  ion  site.  The  diagonalization  involves  the  1 1 
lowest  4/3  \SL\J  manifolds  for  Nd3+  and  the  12  lowest  4/ 7  | SL \J  mani¬ 
folds  for  Gd3+.  The  assumption  is  made  that  Nd3+  substitutes  into  Gd3+ 
sites  of  D2  point  symmetry  in  the  garnet  lattice.  The  CEF-splitting  calcu¬ 
lation  for  trivalent  chromium  assumes  that  the  ions  substitute  into  scan¬ 
dium  (C3|)  sites.  Analyses  of  fluorescence  data  [17]  and  vibronic  spectra 


are  used  to  establish  CEF  levels  of  Cr3+  since  zero-phonon  electronic 
electric-dipole  transitions  are  forbidden  between  3d7  states  at  ion  sites 
having  inversion  symmetry.  The  complete  diagonalization  of  the  3 d 
states  in  a  strong  crystalline  electric  field  of  C3l  symmetry  with  spin-orbit 
interaction  yields  energy  levels  and  wavefunctions  different  from  those 
obtained  from  the  Sugano  and  Tanabe  calculations  for  cubic  symmetry 
with  no  spin-orbit  corrections  [18,19]. 


2.  Experimental  Results 


Single  crystals  were  grown  parallel  to  the  (111)  direction  by  the  Czo- 
chralski  technique.  The  starting  chemicals  were  all  99.999  percent  pure. 
An  indium  crucible  3  in.  in  diameter  and  3  in.  deep  was  used.  The  crystal 
was  rotated  at  15  rpm  and  pulled  at  0.015  in./hr.  The  growth  was  carried 
out  under  an  atmosphere  of  N2  containing  2-percent  02  by  volume;  0-,  is 
added  in  order  to  minimize  the  evaporation  of  gallium  during  growth. 
The  resulting  boule  was  1.5  in.  in  diameter  and  6  in.  long.  A  section  for 
optical  measurement  was  cut  perpendicular  to  the  growth  axis  from  the 
center  portion  of  the  boule,  where  the  calculated  concentrations  were  1.56 
at.  wt.  %  neodymium  based  on  gadolinium  and  1.2  at.  wt.  %  chromium 
based  on  scandium. 

GSGG  is  a  garnet  derivative  of  gadolinium  garnet  (GGG)  where  oc- 
tahedrally  coordinated  gallium  is  replaced  with  scandium.  The  critical 
growth  behavior  of  GSGG  very  closely  resembles  that  of  GGG.  There 
are,  however,  some  minute  differences  as  may  be  expected.  The  con¬ 
gruent  GGG  composition,  with  which  the  crystals  grow,  contains  1  at.  wt. 
%  gadolinium  ions  in  octahedral  sites  (based  on  octahedral  gallium)  as  is 
evidenced  by  comparison  of  lattice  parameters-a  (GGG)  =  12.376  A  and 
a  (congruent  GGG)  =  12.384  A.*  Such  ’’spillover"  of  gadolinium  into  oc¬ 
tahedral  sites  is  not  the  case  in  GSGG.  The  doping  of  GSGG  with 
neodymium  and  chromium  is  better  controlled  as  compared  with 
Y3A1sO)2  (YAG),  The  dopant  distribution  is  governed  by  the  same  prin¬ 
ciples,  but  the  distribution  coefficients  are  much  more  favorable  to  uni¬ 
form  doping.  The  neodymium  distribution  coefficient  K  =  0.6,  and  the 
chromium  distribution  coefficient  approaches  unity  very  closely.  This  re¬ 
sults  in  only  a  very  small  variation  of  neodymium  concentration  in  the 
crystal,  and  the  chromium  distribution  is  uniform. 

The  Ca2+  impurity  concentration  was  below  3  ppm  by  weight  of  scandium 
oxide.  This  is  very  important  with  respect  to  the  absorption  spectrum  of 
Cr3+  in  garnets,  because  the  presence  of  divalent  calcium  affects  the  ab¬ 
sorption  of  chromium  in  garnets  by  changing  relative  intensities  of  the  4A2 
— »  4T2  and  4A2  — >  *TX  transitions,  as  well  as  introducing  a  broad  absorp¬ 
tion  feature  in  the  region  of  0.90  to  1.2  jj.m.  The  nature  of  this  effect  is 


*(A)  +  10  *  ( nm ). 
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not  yet  explained.  The  crystal  grown  for  this  work  contained  such  small 
Ca2+  concentration  that  its  effect  may  be  neglected. 

Absorption  spectra  recorded  between  6000  and  1500  cm  1  with  a  resolu¬ 
tion  of  0.5  cm  were  obtained  from  a  Nicolet  Model  7199  Fourier  trans¬ 
form  infrared  (FTIR)  spectrometer.  Data  were  also  recorded  between 
35, (XX)  and  4,000  cm-1  using  a  Cary  Model  17D  spectrophotometer.  In 
the  ultraviolet  region  the  accuracy  of  the  instrument  is  about  3  A  and  the 
resolution  is  better  than  2  A.  Some  examples  of  this  resolution  are  found 
in  the  figures  representing  the  observed  splitting  of  the  [fiPy  manifolds  of 
Gd3+  in  the  ultraviolet.  The  precision  in  measuring  the  separation  be¬ 
tween  peaks  associated  with  a  given  manifold  at  a  given  temperature  is 
better  than  0.2  A.  In  the  region  of  overlap  between  the  spectrophotometer 
and  the  FTIR  spectrometer,  the  agreement  is  within  one  wavenumber. 

A  conduction  dewar  was  used  to  obtain  spectra  at  liquid  nitrogen  and  liq¬ 
uid  helium  temperatures.  Sample  temperatures  were  not  measured;  the 
spectra  were  recorded  no  sooner  than  30  minutes  after  the  dewar  was 
filled  to  allow  for  thermal  equilibration. 
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3.  Results  and  Conclusions 


Observed  Spectra 

Below  1700  cm  1  the  garnet  lattice  absorbs  [14-16].  Between  1700  and 
6000  cm  the  only  absorption  found  is  associated  with  the  [4/|  1 1/2, 
1 4/ 1 1 3/2,  and  14/]  1 5/2  manifolds  of  Nd3+.  Temperature-dependent  spectra 
were  used  to  establish  the  excited  Stark  levels  of  the  ground  state  mani¬ 
fold  [4/|9/2  of  Nd3+  as  follows:  Z2  =  106  cm  l,Z3  =  167  cm-1,  and  Z4  = 
264  cm  .  Table  1  lists  the  experimentally  established  energy  levels  in  the 
near-infrared  region  obtained  from  an  analysis  of  spectra  recorded  at  three 
different  temperatures.  The  energy  levels  are  similar,  with  few  excep¬ 
tions,  to  those  derived  for  fluorescence  data  by  Kaminskii  and  his 
coworkers  for  Nd3+:GSGG  at  liquid  nitrogen  temperature  [  1 1. 

The  absorption  spectrum  of  the  [4/]  1 1/2  and  [4/]  1 3/2  manifolds  of  Nd3+  at 
liquid  helium  temperature  is  shown  in  figures  1  and  2.  A  sharp  relatively 
weak  peak  appears  to  the  low-energy  side  of  peaks  K,,  Y2,  T3,  Xv  and  Xv 
Peak  Xd  in  figure  2  has  an  unresolved  shoulder  on  each  side.  The  weak 

^  .  —I 

peaks  are  found  within  3  to  8  cm  of  the  strong  absorption  peaks.  With 
the  relatively  high  concentrations  of  Nd3+  and  Cr*+  it  is  not  surprising  that 
some  peak  broadening  and  additional  structure  may  be  observed  under 
high  resolution  fI4|.  Struve  and  Huber  |!7|  interpret  line  broadening  in 
Cr3+:GSGG  as  due  to  lattice  sites  occupied  by  ions  of  different  ionic  radii. 
The  strongest  transitions  appearing  in  figures  1  and  2  have  linewidths  at 
half-maximum  less  than  10  cm  1  at  liquid  helium  temperature. 

There  is  no  absorption  above  the  1 4/ 1 1 5/2  manifold  until  one  reaches  the 
manifolds  of  14FJ3/2  (1 1,450  cm"1)  (fig.  3),  [4FJ5/2  and  [2H\9/2  (12,500 
cm  '),  and  [4F]7/2  and  |4S]3/2  (13,500  cm  *)•  The  excited  Stark  levels 
for  each  of  these  manifolds  have  been  identified  from  the  liquid  helium 
temperature  spectrum  (table  2).  Temperature-dependent  spectra  were 
used  to  establish  the  ground  state  splitting  of  the  1 4/J9/2  manifold  consis¬ 
tent  (within  a  wavenumber)  to  the  ground  state  splitting  established  from 
the  temperature-dependent  spectra  of  the  [4/]  11/2,  [4/]  13/2,  and  14/]  1 5/2 
manifolds  (table  1). 

Absorption  by  Cr3+  appears  near  7000  A  (14,200  cm  ').  Only  at  liquid 
helium  temperature  are  the  details  sharp  enough  to  measure  (table  3).  The 
/?,  and  R2  peaks  of  the  |2£]  state  (levels  3  and  4  in  fig.  4)  are  found  at 
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6965  A  (14,354  cm  ')  and  6951  A  (14,382  cm  *).  Both  peaks  have  a 
linewidth  of  10  A  at  half-maximum.  R2  (level  4)  shows  evidence  of  in¬ 
homogeneous  line  broadening  (fig.  4).  A  band  partially  resolved  into 
three  peaks  nearly  as  intense  as  the  /?,  and  R2  peaks  lies  to  the  high- 
energy  side  of  the  R2  peak.  The  peaks  are  at  6933  A  (14,420  cm  '),  6923 
A  (14,440  cm  1 ),  and  6909  A  (14,470  cm  *).  The  onset  of  the  first  strong 
absorption  band  of  Cr  begins  around  6900  A  and  shows  a  distinct  shoul¬ 
der  near  6632  A.  This  shoulder  has  structure  with  peaks  at  6879,  6795, 
and  6776  A  (table  3).  The  band  has  a  maximum  absorbance  at  6360  A 
which  decreases  to  the  band  edge  at  5500  A.  Superimposed  on  the  first 
Cr3*  band  are  the  absorption  peaks  of  the  Nd3*  manifolds  |4F|9/2  (fig.  4), 
\2H\  1 1/2  (fig.  5),  (2G|7/2,  and  |4G|5/2.  Analysis  of  the  absorption  spec¬ 
trum  of  Cr3+:GSGG  (without  Nd3*)  and  GSGG  (without  Cr3*  and  Nd3*) 
allows  one  to  separate  spectral  details  of  Cr3+  and  Nd3+.  The  absorption 
spectrum  of  Cr3+:GSGG  is  displayed  in  figure  6. 

Nd3*  absorption  associated  with  the  [4G|7/2,  [2GJ9/2,  and  |2Af|13/2  mani¬ 
folds  appears  near  the  onset  of  the  second  strong  Cr3*  band  beginning  at 
5150  A.  The  Nd3*  energy  levels  for  these  states  are  listed  in  table  2.  Near 
the  beginning  of  the  second  band,  Cr3*  peaks  appear  at  5043  and  4993  A 
(table  3).  These  peaks  have  linewidths  and  line  shapes  similar  to  the  R2 
and  /?,  peaks.  The  second  Cr3*  band  is  composed  of  two  peaks  at  4806 
and  4540  A.  The  band  edge  on  the  high-energy  side  is  at  3900  A.  Super¬ 
imposed  on  the  second  Cr3*  band  are  the  Nd3*  spectra  for  the  |4GJ9/2, 
|4G|1 1/2,  [2ATJ 1 5/2,  [2P]l/2,  and  1 2D] 5/2  manifolds.  Table  2  presents  the 
remaining  energy  levels  for  Nd3+,  including  those  found  on  a  third  Cr3* 
band  beginning  at  3400  A  (fig.  7).  Cr3*  peaks  identified  from  Cr3*:GSGG 
(no  Nd3*)  include  3665,  3438,  3431,  3425,  and  3395  A.  The  peaks  have 
been  confirmed  by  scanning  the  spectrum  of  the  undoped  crystal.  Cr3* 
absorption  peaks  on  the  third  Cr3*  band  are  found  at  3198,  3119,  3116, 
3105,  3100,  3085,  and  3081  A  (table  3).  These  peaks  are  weak  in  com¬ 
parison  with  the  strong  sharp  spectra  of  Gd3*  (fig.  8).  Below  2800  A  the 
absorption  of  the  crystal  is  too  intense  to  permit  observation  of  additional 
structure.  Table  4  lists  the  energy  levels  of  Gd3*  measured  at  three  differ¬ 
ent  temperatures.  Even  at  room  temperature  it  is  possible  to  see  all  ex¬ 
pected  Stark  levels  resolved  in  each  \6P\J  manifold  (fig.  8). 


Table  I.  Infrared  energy  levels:  Nd3* 


Value  of  £  (cm-1) 
at  various  temperatures0 


Label  RT 


Calculated 

£  A£ 

(cm-V  (cm-1)c 


Spectra  recorded  using  the  FTIR  spectrometer  at  room  temperature  (RT),  liquid 
nitrogen  (LN)  temperature,  and  liquid  helium  (LHe)  temperature.  141 1 1512 
manifold  recorded  on  spectrophotometer.  Where  data  overlap,  energy  levels 
agree  within  /  cm  . 

bCalculated  energy  levels  based  on  parameters  appearing  in  table  6  under 
column  6.  Nd* '(B). 

c Difference  between  calculated  and  observed  levels  at  Lllc  temperatures. 
dFluorescencc  spectra  at  LN  temperature  //  /  provide  complete  splitting  of 
141 1912  manifold  as  follows:  Z;  »  0.  Z2  »  107,  Z}  «  168.  Vu  *  263,  Z3  »  763,  all 


‘ Theoretical  centroid  in  cm~ 


Figure  1.  Absorption 
spectrum  of  |4/|l  1/2 
manifold  of  Nd3*  recorded 
at  liquid  helium  tempera¬ 
ture  with  FTIR 
spectrometer. 


Figure  2.  Absorption 
spectrum  of  |4/]  1 3/2 
manifold  of  Nd3*  recorded 
at  liquid  helium  tempera¬ 
ture  with  FTIR 
spectrometer. 


Figure  3.  Absorption 
spectrum  of  |4F|3/2  mani¬ 
fold  of  Nd3*  recorded  at 
room  temperature  showing 
temperature-dependent 
transitions  from  Stark  lev¬ 
els  Zj  (106  cm"1),  Zj  (167 
cm  '),  and  Z4  (263  cm'1) 
of  ground  state  manifold 
I4/J9/2. 


\SL\J 

Label 

Value  of  £  (cm  ) 
at  various  temperatures'1 

RT  LN  LHe 

Calculated 

£ 

(cm-1)* 

m\m 

11,431 

11,432 

11,432 

11,432 

11, 497* 

/?2 

11,492 

11,494 

11,494 

11,494 

[4f]  5/2 

Si 

12,350 

12,352 

12,354 

12,348 

12,4 74* 

Si 

12,395 

12,390 

12,390 

12,392 

12// 19/2 

Si 

12,395 

12,404 

12,405 

12,405 

1 2,592-' 

S 4 

12,540 

12,543 

12,544 

12,544 

s5 

12,580 

12,584 

12,586 

12,590 

*6 

12,610 

12,617 

12,618 

12,618 

Si 

12,743 

12,745 

12,745 

12,746 

s$ 

12,790 

12,791 

12,793 

12,791 

mm 

*1 

13,372 

13374 

13,374 

13,375 

13.455“' 

a2 

13,429 

13,430 

13,430 

13,432 

13,550 

13353 

13355 

13356 

(4513/2 

13,560 

13,566 

13,566 

13367 

13.547“' 

As 

13,560 

13,572 

13.573 

13,572 

*6 

13,608 

13,610 

13,611 

13,609 

(4F19/2 

Bi 

14.634 

14,637 

14,639 

14,641 

14,724“' 

b2 

14,669 

14,671 

14,673 

14,675 

*3 

14,786 

14,789 

14,791 

14,789 

B 4 

14,804 

14,806 

14,810 

14,809 

*5 

14,899 

14,900 

14,902 

14,901 

[2//JU/2 

C, 

15,864 

15,865 

15,865 

C2 

15,888 

15,890 

15,890 

C3 

15,965 

15,966 

15,967 

c. 

16,086 

16,090 

16,091 

16,105 

16,109 

16,110 

(4CJ5/2 

£>. 

16,907 

16,909 

16,910 

02 

17,000 

17,005 

17,008 

12CI7/2 

£>3 

... 

17,020 

17,022 

£>< 

17,070 

17,074 

17,075 

*>5 

17,278 

17,282 

17,284 

£>6 

... 

17350 

17357 

d7 

17,563 

17370 

17373 

Table  2.  Visible  and  ultraviolet  energy  levels:  Nd?+  (cont’d) 


Value  of  £  (cm  ) 
at  various  temperatures4 


Value  of  £  (cm  ) 
at  various  temperatures4 


l  SLU 


I4CI7/2 


|2/V|  1 3/2 


I2CI9/2 


[4GJ9/2 


[4GJI1/2 


[2/T1 1 5/2 


Label 

RT 

LN 

LHe 

[SLU 

Label 

RT 

LN 

LHe 

18,777 

18,781 

18,783 

|2/f  115/2 

"s 

21,778 

21,780 

21,783 

18,864 

18,864 

18,864 

(cont’d) 

21,884 

21,884 

21,885 

*3 

18,881 

18,883 

18,884 

"7 

22,004 

22,005 

22,006 

18,997 

18,901 

19,004 

22,053 

22,054 

22,054 

19,203 

19,205 

19,207 

|2/»|l/2 

A 

23,190 

23,192 

23,193 

19,236 

19,238 

19,238 

[2£»]5/2 

23,723 

23,723 

23,724 

19,330 

19,334 

19,335 

'2 

23,778 

23,779 

3,780 

F* 

19,354 

19,359 

19,359 

•A 

23,832 

23,840 

23,843 

F> 

19,410 

19,411 

19,412 

f6 

19,480 

19,487 

19,488 

[2/>J3/2 

26,045 

26,046 

26,047 

Fj 

19,518 

19,520 

19,521 

*2 

... 

26,055 

26,056 

Ft 

19,586 

19,590 

19,592 

f9 

19,615 

19,619 

19,620 

[4D]3/2 

27,685 

27,689 

27,689 

F>  0 

19,670 

19,672 

19,672 

^2 

27,761 

27,765 

27,765 

^1. 

19,820 

19,824 

19,824 

F\2 

19,890 

19,895 

19,895 

[4015/2 

27,880 

27,880 

27,882 

28,035 

28,039 

28,040 

Gt 

20,793 

20,795 

20,797 

c,i 

20,824 

20,828 

20,828 

[4D11/2 

28,105 

28,112 

28.113 

c3 

20,842 

20,845 

20,846 

Af4 

28,145 

28,150 

28,159 

Afs 

28,280 

28,283 

28,284 

20,852 

20,856 

20,858 

Gs 

20,977 

20,980 

20,980 

[2/J11/2 

*6 

28328 

28333 

28335 

21,040 

21,042 

21,043 

Af7 

28,409 

28.413 

28,413 

G7 

21,093 

21,096 

21,096 

M, 

28,705 

28,712 

28,714 

C, 

21,143 

21,145 

21,146 

Af, 

28,853 

28.857 

28.859 

G, 

21,158 

21,159 

21,159 

*.0 

29,014 

29,017 

29,018 

^10 

21,203 

21,206 

21,207 

c„ 

21,234 

21,238 

21,240 

(2LJ15/2 

A'. 

29,853 

29,854 

29,856 

A'j 

29,963 

29,960 

29,958 

"■ 

21,609 

21,610 

21,611 

A', 

... 

... 

29,969 

21,714 

21,714 

21,714 

"3 

21,740 

21,740 

21,742 

[40)7/2 

29,993 

29,997 

30,000 

"4 

21,772 

21,773 

21,774 

30,051 

30,056 

30,059 

V 

s’ 


Table  2.  Visible  and  ultraviolet  energy  levels:  Nd3*  (conl'd) 

Value  of  E  (cm"1)  Value  of  £  (cm  ') 

at  various  temperatures*1  at  various  temperatures** 


1 SLU 

Label 

RT 

LN 

LHc 

[SL]J 

Label 

RT 

LN 

LHe 

140)7/2 

"• 

30,070 

30.073 

[2LJ17/2 

o. 

31,413 

31,413 

31,414 

(conl’d) 

"7 

30,178 

30,180 

30,180 

02 

31,518 

31,520 

31,522 

", 

30,209 

30,210 

30,210 

Os 

31,568 

31,570 

31,570 

", 

30,263 

30,264 

30,265 

o4 

31,610 

31,610 

31,612 

",o 

30,300 

30,306 

30,308 

Os 

31,660 

31,662 

31,662 

30,352 

30,356 

30,359 

o6 

31,751 

31,753 

31,755 

12/113/2 

",  2 

30,428 

30,430 

30,431 

[2//J9/2 

Py 

32,633 

32,633 

32,634 

30,453 

30,453 

30,454 

Pi 

32,684 

32,686 

32,686 

"m 

30,506 

30,509 

30,510 

P3 

32,708 

32,710 

32,713 

"|5 

30,532 

30,536 

30,540 

Pa 

32,770 

32,772 

32,774 

",6 

30.625 

30,627 

30,630 

PS 

32,791 

32,795 

32,797 

“Spectra  recorded  on  spectrophotometer  at  room  temperature  ( RT ),  liquid  nitrogen  (LN)  temperature,  and  liquid 
helium  (Llle)  temperature. 

bCalculatcd  energy  levels  based  on  parameters  appearing  in  table  6  under  column  6,  Nd’^fB). 
c Difference  between  calculated  and  observed  levels  at  Llle  temperature;  rms  deviation  for  45  levels  (tables  I  and  2) 
is  1.7  cm'1. 


Bandwidth:  full  width  at  half-maximum. 

bThe  parameters  (in  cm  ')  used  in  the  calculation  are  Fl2>  =  54320,  F*4)  «  45.094,  a  -  2.88,  y  « 
-63.28,  ^  -  169.64,  B20  =  1072,  B40  =  -22251,  and  B4}  =  23,443,  with  an  rms  =  87.7  cm'\  ( B  = 
620.  C  -  3420.  Dq  >  /450.  v  =  -577.5,  v '  =  5123). 

cThe  irreducible  representations  are  4  =  V4  +  and  6  =  V6  (doublet)  for  the  group  C}  or  C}i  ( all  r 
and  T, ).  Ihe  notation  is  from  Koster  et  al  I30j. 
d()nly  states  of  I  percent  or  greater  are  listed. 

* The  splitting  of  the  ground  stale  is  0.235  cm~]  with  small  amounts  of  the  *P  and  2G  mixed  in. 

1 Shoulder  of  a  band. 


Figure  4.  Absorption 
spectrum  of  Cr3*  levels  3 
through  10  and  the 
I4FI9/2  manifold  (B 
group)  of  Nd3*  recorded  at 
liquid  helium  temperature. 


Figure?.  Absorption 
spectrum  of  (2//IU/2 
manifold  of  Nd3*  superim¬ 
posed  on  first  strong  Cr3* 
band,  recorded  at  liquid 
helium  temperature. 
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Figure  n.  Absorption 
spectrum  of  Cr^rGSGC 
between  3000  and  8000  A 
recorded  at  room  tempera¬ 
ture  with  a  Cary  Model  17 
spectrophotometer. 


WAVELENGTH.  A 


Figure  7.  Absorption 
spectrum  of  (40 1 3/2  (L 
group),  (40(5/2  and 
(2/(1 1/2  (M  group),  and 
part  of  (40(7/2  and 
(2/|13/2  (N  group)  mani¬ 
folds  of  Nd3*  superim¬ 
posed  on  third  Cr3*  band 
recorded  at  liquid  helium 
temperature. 
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Figure  8.  Absorption 
spectrum  of  [6P|7/2|  (A 
group),  |6P|5/2  (fl  group), 
and  |6/*|3/2  (C  group) 
manifolds  of  Gd3*  re¬ 
corded  at  room 
temperature. 


WAVELENGTH,  A 


Table  4.  Ultraviolet  energy  levels:  (Jd 


“Spectra  recorded  on  spectrophotometer  at  room  temperature  (RT), 
liquid  nitrogen  ILN)  temperature,  and  liquid  helium  (Lite)  temperature. 
kCalc.ulaied  energy  levels  bated  on  parameters  appearing  in  table  6. 
c Difference  between  calculated  and  observed  levels  at  Lllc  temperature; 
rms  deviation  for  14  levels  is  4  cm~' . 
d Theoretical  centroid  in  cm\ 


3.2  Crystal-Field  Splitting  Calculations:  Nd3+,  Gd3+ 

The  material  GSGG  is  a  mixed  garnet  with  the  same  crystal  structure  as 
Y3A150,2  (YAG)  (21,22).  The  cubic  space  group  is  laM  with  site  sym¬ 
metries  as  follows:  Gd(Nd),  D2;  Sc(Cr),  C3(;  Ga,  S4;  and  O,  Cj.  Crystal¬ 
lographic  data  are  summarized  in  table  5. 

The  CEF  splitting  of  the  4/"  [SL|7  manifolds,  assuming  that  Gd3+  and 
Nd3+  occupy  sites  of  D2  point  symmetry,  was  analyzed  by  diagonalizing  a 
parameterized  Hamiltonian, 


Hl)i  =  2.  Bbn  Clan  ,  (1) 

km 

in  a  free-ion  wave  function  basis  involving  the  1 1  lowest  4 f  | SL\.f  mani¬ 
folds  for  Nd3+  and  the  12  lowest  4/7  |SL|7  manifolds  for  Gd3+  (23|. 
Free-ion  wave  functions  and  reduced  matrix  elements  of  the  l/k)  unit 
spherical  tensors  were  first  obtained  by  diagonalization  of  the  free-ion 
Hamiltonian  using  parameters  given  by  Carnall,  Fields,  and  Rajnak  |24]. 
The  nine  real,  even-fold  (even-*)  parameters  in  the  D2  Hamiltonian  were 
varied  to  obtain  agreement  with  observed  splitting.  The  relationship  be¬ 
tween  the  8^  parameters  used  in  this  study,  along  with  the  A”  (r* )  and 
Bkm  parameters  which  are  frequently  reported,  and  the  six  possible  equiva¬ 
lent  sets  of  Bfa  parameters  for  D2  symmetry  can  be  obtained  by  using  the 
relations  given  elsewhere  [14J. 


Table  5.  Crystallographic  and  x-ray  data  of  Gd3Sc2Ga30I2* 
(Cubic  laid,  230,  Z  =  8) 


Ion 

Site 

Symmetry 

X 

y 

z 

Gd 

24c 

0 

1/4 

1/8 

Sc 

16a 

c* 

0 

0 

0 

Ga 

24  d 

0 

1/4 

3/8 

O 

96h 

c, 

-0.0272 

0.0558 

0.1501 

* Lattice  constant  is  a  =  12.5668  A  from  Brandle  and  Barns  138 / 
and  the  fractional  positions  for  the  oxygen  ions  are  for  Y}Ga}On 
from  Euler  and  Bruce  I39j. 


1 


I 


ri 


a 

‘V 
’  V 
1*1 1 
Vt 


a 


a 


a 


An  initial  set  of  phenomenological  parameters  for  Nd3+  was  obtained 
from  Morrison,  Leavitt,  and  Gildner  |23),  who  calculated  the  CEF  split¬ 
ting  of  the  14/J./  and  [4FJ3/2  manifolds  using  the  data  from  Kaminskii  et 
al  1 1 J.  The  present  paper  extends  that  work  by  including  data  and  calcula¬ 
tions  for  the  additional  Nd3+  manifolds,  [4FJ5/2,  [2WJ9/2,  |4F]7/2, 
I45J3/2,  and  \4F\9J2.  Table  6  presents  the  Final  set  of  parameters 
based  on  the  experimental  energy  levels  listed  in  tables  I  and  2.  The  cal¬ 


culated  splittings  are  also  listed  in  tables  1  and  2.  An  mis  deviation  of 
1 .7  cm  is  obtained  for  45  observed  levels  using  the  Final  set  of 


phenomenological  B ^  parameters  under  column  Nd  (B)  (table  6). 


The  phenomenological  B ^  parameters  were  compared  with  those  ob¬ 
tained  from  a  lattice  sum  calculation  by  Morrison,  Leavitt,  and  Gildner 
[23],  The  fractional  oxygen  positions  were  taken  from  the  isostructural 
material  Y3Ga5012  |14|.  Effective  charges,  Z-  (qi  =  eZ)  in  units  of  the 
magnitude  of  the  electric  charge,  were  taken  as  Z0d3-  =  Z^.  =  3,  Z0,j-  =  1 , 
and  Z02  =  -1.5.  The  polarizability  of  oxygen,  o^,  was  taken  as  0.244 
A3.  The  resulting  lattice  sums,  A ^  which  include  the  point-charge, 
point-dipole,  and  self-induced  contributions,  are  listed  in  table  6.  The  pa¬ 
rameters  A^  are  related  to  the  phenomenological  parameters  through 
the  expression 


Bkm  =  p*  (Nd)  A  km  » 


where  pt(Nd)  are  radial  factors  given  by  Morrison,  Karayianis,  and 
Wortman  125).  The  lattice  sum  calculation  predicts  reasonable  values  for 
A  km  considering  the  assumptions  that  have  been  made. 


The  set  of  B ^  parameters  for  Gd3+  was  obtained  by  using  pt(Gd)  [25 1  and 
the  phenomenological  A ^  lattice  parameters  from  table  6.  The  number 


(14)  of  observed  energy  levels  for  Gd  is  marginal  for  a  meaningful 


analysis,  when  the  symmetry  is  low  and  requires  nine  CEF  parameters. 
The  rms  deviation  for  14  observed  levels  is  4  cm-1,  which  is  within  the 
uncertainty  associated  with  the  data  recorded  in  the  ultraviolet  region. 


V.**' 


& 

Table  6.  CEF  parameters  for  Nd',+  and  Gd3 


Alcaic.)"  A^Cphenom.)* 

A.  (cm'/A)"  Nd  (cm’ ’/A)"  B. 


Nd3+(Af  Nd3+(B)rf 
(cm-1)  (cm'1) 


(cm'1) 


“Table  2,  ref.  23. 

b\'alues  of  parameters  calculated  from  phenomenological  parameters  given  in 

ref.  23.  an  analysis  of  fluorescence  data  for  [4IJJ  and  [4113/2  manifolds  reported  in  ref.  I. 
initial  set  of  parameters,  from  ref.  23. 

dFinal  set  of  B^  parameters  for  Nd 3t ;  45  levels;  rms  deviation  / .  7  cm'1 . 

‘B^  parameters  for  Gd 3t ;  14  levels;  rms  deviation  4  cm. 

3.3  Crystal-Field  Splitting  Calculations:  Cr3+ 

Since  we  have  not  previously  given  a  detailed  description  of  the  method 
used  to  calculate  the  energy  level  and  the  fitting  procedures  used  in  the 
theory  of  transition-metal  ions  in  solids,  we  go  into  some  detail  here.  The 
free-ion  interactions  ( Hri )  we  include  are  as  follows: 

N  N 

Hr,  =  I  /**>  lC^(OC*0')  +  aL(£+l)  +  YG(/?5)  +C*  I  /,  ■  s.  ,  (3) 

*=2,4  i>j  i«l 

where  F*k)  are  the  Slater  integrals,  a  and  y  are  parameters  of  the  Trees 
interactions,  and  is  the  spin-orbit  constant.  In  equation  (3),  the  Slater 
integrals  are  related  to  the  Racah  parameters  by 


F<2)  =  1(1B  +  C)  , 
f*4)  =  63C/5  , 

and  the  are  related  to  the  spherical  harmonics  by 

C*  (/)  =  Mn/(2k  +  1)  Ytu,  (0, ,  <(>0  , 

Ck-<,  =  (-!)’  C;  . 


The  quantity  G(/?5)  is  the  Casimir  operator  for  the  rotation  group  Z?5,  and 
the  values  for  all  the  states  of  dN  are  given  by  Judd  |26J.  For  the  Cr3+ 
free  ion,  the  parameters  in  equation  (3)  have  the  following  values  from 
Uylings  et  al  [27 1:  F<2)  =  74,201,  A7*4’  =  45,822,  a  =  29.87,  and  £ \d  =  275 
(all  in  units  of  cm  '). 

Using  equation  (3)  in  a  least-squares  fit  to  the  data  reported  by  Sugar  and 
Corliss  [281,  we  obtained  the  following  values:  F*2’  =  72,389,  F(4)  = 
43,044,  a  =  91.75,  y  =  -129.6,  and  ^  =  277.7  (all  in  units  of  cm  '). 

The  difference  in  the  two  sets  of  parameters  is  mostly  because  Uylings  et 
al  I27j  include  the  Trees  interaction  with  the  parameter  T  as  well  as  P, 
which  multiplies  the  seniority  operator.  Because  of  the  agreement  given 
above  for  transition-metal  ions  in  solids,  equation  (3)  can  be  considered  an 
adequate  representation  of  the  free-ion  interaction. 

For  threefold  symmetry  such  as  is  found  at  the  Sc  site  in  GSGG  (see  table 
6),  the  crystal-field  interaction  ( HCF )  is  taken  as 

N  N  N 

Hcf  -  B20  Z  C20  (1)  +  Bm)  Z  C40  (i)  +B43  Z  [C4-3O)  -  C43  (()]  ,  (6) 

i»l  i=»l  1=1 

where  the  B ^  are  the  crystal-field  parameters  and  fi43  can  be  chosen  real 
and  positive  with  no  loss  in  generality.  In  the  cubic  approximation,  fljn  = 
0  and  fl43  =  V 10/7  |A?4o|,  and  for  this  symmetry  (threefold),  Bm  =  -14 Dq, 
where  we  assume  Dq  is  a  positive  number. 

The  matrix  elements  of  the  Hamiltonian  given  in  equations  (3)  and  (6) 
were  computed  using  tapes  of  coefficients  of  fractional  parentage  for  the 
(P  configuration.  The  notation  of  the  free-ion  levels  is  that  of  Nielson  and 
Koster  129].  Total  angular  momentum  wave  functions  are  used;  for  C3 
symmetry  the  T4  matrix  elements  [30]  are  obtained  by  choosing  the 
projection  My,  as  Mj=  1/2  +  3 q,  with  q  any  integer  and  I  Mj\  <J.  For  the 
r6  irreducible  representation,  again  we  have  Mj  =  3/2  +  3 q,  where  q  is 
any  integer  and  \Mji  <  J.  With  these  values  of  Mj  the  basis  is  chosen  as 
Mj )  +  (—  1  \J—M j )  I'll  and  the  energy  levels  for  the  basis 
\JMj)  -  I'll,  are  identical.  With  these  basis  functions. 


the  matrices  for  the  entire  Hamiltonian  are  39  x  39  for  T4  (r5  is 
degenerate  with  r4)  and  21  x  21  for  T6.  No  a  priori  assumptions  are 
made  concerning  the  relative  strengths  of  the  various  terms  in  the  Hamil¬ 
tonian.  This  choice  of  basis,  while  not  efficient  for  a  cubic  environment, 
avoids  the  tedious  perturbation  calculations  when  the  ion  is  definitely  not 
in  a  site  of  cubic  symmetry. 

Having  established  the  method  by  which  the  calculations  are  made,  we 
now  turn  briefly  to  describe  the  method  used  to  establish  the  experimental 
energy  levels.  Since  C3l  symmetry  has  inversion,  pure  electronic 
electric-dipole  (zero-phonon)  transitions  are  forbidden  between  states  of 
the  same  parity.  However,  because  the  coupling  between  the  lattice  and 
the  electronic  states  is  so  strong,  the  vibronic  spectra  are  intense  and  rep¬ 
resent  most  of  the  observed  Cr3+  absorption  spectra.  The  only  exception 
is  the  observation  of  several  allowed  electronic  magnetic-dipole  transi¬ 
tions  between  Stark  levels  that  are  nearly  100-percent  pure  V  suites  such 
as  transitions  between  the  ground  state  level  and  levels  3  and  4  (table  3). 

At  liquid  helium  temperature  most  of  the  temperature-dependent  vibronic 
spectrum  is  diminished,  so  that  if  a  zero-phonon  or  electronic  transition 
were  allowed,  it  would  be  found  near  the  low-energy  side  of  the  observed 
vibronic  spectra  associated  with  an  excited  Stark  level.  Some  years  ago 
Satten  and  his  coworkers  were  successful  in  analyzing  the  vibronic 
spectra,  especially  of  octahedrally  coordinated  uranium  complexes  1 3 1 1. 
The  method  of  analysis  they  used  was  adopted  and  expanded  by  others 
analyzing  vibronic  spectra  of  rare-earth  ions  having  different  symmetries 
[32|.  Recently  vibronic  states  in  C3l  symmetry  were  analyzed  and  the 
electronic  energy  levels  reported  |33).  The  method  of  analysis  of  vibronic 
spectra  used  here  follows  the  detailed  descriptions  reported  earlier.  We 
conclude  from  our  assignments  that  the  uncertainty  in  location  of  isolated 
excited  Stark  levels  is  within  the  final  rms  value  reported.  The  fluores¬ 
cence  spectrum  observed  at  liquid  helium  temperature  [171  also  helps  to 
establish  the  location  of  the  levels.  Several  authors  described  the  use  of 
fluorescence  data  to  establish  electronic  energy  levels  1 14,17,31 1. 

To  determine  approximate  atomic  parameters,  we  began  by  assuming  cu¬ 
bic  symmetry.  From  observed  spectra,  we  estimated  the  position  of  the 
experimental  centroids  as  4A2  =  0,  2E  =  14,368,  ^T,  =  14,436,  = 

14,793,  47' ,  =  20,330,  and  4T^  =  32,251  (in  cm  *).  In  making  these  selec- 


V 

I1 


tions,  we  used  a  Tanabe  plot  [34]  with  Dq  =  1479.3.  This  value  of  Dq 
gave  B40  =  -20,710  cm  1  and  fl43  =  24,753  cm  '.  For  free-ion 
parameters,  we  started  with  those  of  Cr3+  in  Y3A15012  (YAG)  [35],  which 
are  F{1)  =  55,800  cm  ',  F*4)  =  36,800  cm'1,  and  a  =  y  =  =  0.  With  the 

above  parameters  as  starting  values,  a  least-squares  fit  to  the  experimental 
centroids  was  performed.  The  resulting  parameters  obtained  were  F(2>  = 
50,200,  F(4)  =  40,821,  a  =  35.09  (y  =  (,d  =  S20  =  0  not  varied),  fi40  = 
-20,797,  and  B43  =  24,857.71  (in  cm  ').  A  similar  fit  to  the  limited  data 
of  Struve  and  Huber  [17|  gave  F(2)  =  53,961,  F<4)  =  40,782,  Bw  = 
-2 1 ,882,  and  B43  =  26, 1 54  (in  cm  1 )  (a  =  y  =  C,d  =  fl2 o  =  0). 

Proceeding  now  with  the  correct  symmetry  (C3-)  and  the  data  given  in  ta¬ 
ble  5,  we  computed  the  point  charge  lattice  sum  parameters  A ^  ]23], 
which  are  A20  =  812  citTVA2,  A40  =  -1 1,413  cm  '/A4,  and  A43  =  13,327 
cm  /A4.  We  then  used  the  rotational  invariants  136]  defined  here  as 


Sn(B) 


Blo+2  1 


n  -]'« 

£  BnmBnm  , 

»0 


which  for  C*  symmetry  gives 

=  \B>\)  +  2B43],/2  .  (8) 

If  we  assume  the  theoretical  B ^  are  given  by  B^  =  p nAwn,  then  p4  = 
S4(B)/S4(A).  From  equations  (7)  and  (8)  above  we  have 

p4=  1.8537  A4  .  (9) 

If  further  we  assume  that  pn  =  (r")w/r/x  "  (with  (/*)w/7  given  by  a  Hartree- 
Fock  calculation  and  x  a  radial  expansion  parameter),  then  we  can  write 

P2  =  (f2  )HF  Vp 4/(^>HF  •  ^ 

From  Fraga  et  al  [37]  we  have  <r2)///r  =  0.4018  A2  and  (r4);//7  =  0.3344  A4; 
therefore,  p2  =  0.9460  A2.  Using  the  relation  B‘ltm  =  pnA**  and  these  val¬ 
ues  of  p2  and  p4  we  have 

B'20  =  786  cm-1  , 

flio  =-21,156  cm"'  ,  (11) 

=  24,704  cm"'  . 


The  crystal-field  parameters  given  in  equation  (11)  were  then  used  as 
starting  values  along  with  the  free-ion  parameters  obtained  from  fitting  the 
experimental  centroids  and  with  C,d  =  200  cm-1  in  a  least-square  fit  to  all 
the  established  energy  levels.  A  total  of  35  levels  of  Cr3+  were  involved 
in  the  final  analysis.  The  parameters  obtained  were  F42*  =  54,320,  F44*  = 
43,094,  a  =  2.88,  y=  -63.28,  £rf=  169.64,  fl20  =  1,072,  fl40  =  -22,251,  and 
fi43  =  23,443  (in  cm  ).  The  rms  deviation  is  87.7  cm-1. 

Table  3  gives  the  resulting  energy  levels,  along  with  the  percentage  of  the 
free-ion  composition  of  each  state  (only  values  greater  than  1  percent  are 
given).  The  number  of  levels  of  the  free-ion  composition  is  limited  to  3; 
the  labels  are  in  the  convention  of  Nielson  and  Koster  |29].  Because  of 
the  spin-orbit  coupling  and  the  twofold  crystal  field  the  "2£"  cubic- 
level  becomes  predominantly  "4F,"  and  the  entire  region  from  14,300 
through  15,000  cm  is  so  mixed  that  a  cubic  interpretation  is  practically 
impossible. 


From  temperature-dependent  absorption  spectra  of  Nd3+:GSGG  |1], 
Cr3+:GSGG  1 17J,  and  Nd3+:Cr3+:GSGG  it  has  been  possible  to  catalogue 
over  400  zero-phonon  transitions  which  were  used  to  establish  the  Stark 
levels  of  Nd3+  and  Gd3+  ions  occupying  D2  point  symmetry  in  the  lattice. 
The  vibronic  spectrum  of  Cr3+  in  C3(  sites  has  been  analyzed  to  establish 
the  approximate  location  of  many  of  the  Stark  levels  of  Cr+.  A  Hamil¬ 
tonian  consisting  of  Coulombic,  spin-orbit,  and  CEF  terms  and  having 
symmetry  appropriate  to  each  ion  was  diagonalized  to  obtain  theoretical 
energy  levels.  The  rms  deviation  between  calculated  and  observed  levels 
for  Nd3+  (with  45  levels)  is  1.7  cm  ',  for  Gd3+  (with  14  levels)  is  4  cm  1 
and  for  Cr  (with  35  levels)  is  87.7  cm  .  The  lattice  sum  calculations 
for  both  Nd3+  and  Cr3+  are  in  reasonable  agreement  with  the  phenom¬ 
enological  parameters  obtained  from  fitting  the  observed  energy  levels. 
The  C3i  point  symmetry  of  the  lattice  must  be  used  to  interpret  the  energy 
levels  of  Cr3+  in  GSGG. 
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